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ABSTRACT 

We have calculated synchrotron spectra of relativistic blast waves, and find predicted characteristic frequencies 
that are more than an order of magnitude different from previous calculations. For the case of an adiabatically 
Q\ ■ expanding blast wave, which is applicable to observed gamma-ray burst (GRB) afterglows at late times, we give 

Q\ ' expressions to infer the physical properties of the afterglow from the measured spectral features. 

ON . We show that enough data exist for GRB 970508 to compute unambiguously the ambient density, n = 0.03 cm" 3 , 

and the blast wave energy per unit solid angle, £ = 3 x 10 52 erg/47r sr. We also compute the energy density in 
l-H ' electrons and magnetic field. We find that they are 12% and 9%, respectively, of the nucleon energy density and 

Qh, thus confirm for the first time that both are close to but below equipartition. 

For GRB 971214, we discuss the break found in its spectrum by Ramaprakash et al. (1998). It can be interpreted 
either as the peak frequency or as the cooling frequency; both interpretations have some problems, but on balance 
the break is more likely to be the cooling frequency. Even when we assume this, our ignorance of the self- 
absorption frequency and presence or absence of beaming make it impossible to constrain the physical parameters 
^ i of GRB 971214 very well. 

Subject headings: gamma rays: bursts — gamma rays: individual (GRB 970508) — gamma rays: individual 
(GRB 971214) — gamma rays: theory — acceleration of particles — magnetic fields 

m ■ 

' 1. INTRODUCTION 

o ; 

Explosive models of gamma-ray bursts, in which relativistic ejecta radiate away some of their kinetic energy as they are slowed 
. down by swept-up material, naturally lead to a gradual softening of the emission at late times. This late-time softer radiation has 
been dubbed the 'afterglow' of the burst, and its strength and time dependence were predicted theoretically (Meszaros and Rees, 
1997). Soon after this prediction, the accurate location of GRB 970228 by the BeppoSAX satellite's Wide Field Cameras (Piro et al. 
1995, Jager et al. 1995) enabled the detection of the first X-ray and optical afterglow (Costa et al. 1997, Van Paradijs et al. 1997). Its 
behavior agreed well with the simple predictions (Wijers et al. 1997, Waxman 1997a, Reichart 1997). 

The basic model is a point explosion with an energy of order 10 52 erg, which expands with high Lorentz factor into its surroundings. 

■ As the mass swept up by the explosion begins to be significant, it converts its kinetic energy to heat in a strong shock. The hot, 
shocked matter acquires embedded magnetic fields and accelerated electrons, which then produce the radiation we see via synchrotron 
emission. The phenomenon is thus very much the relativistic analogue of supernova remnant evolution, played out apparently in 
seconds due to the strong time contractions resulting from the high Lorentz factors involved. Naturally, the Lorentz factor of the 
blast wave decreases as more matter is swept up, and consequently the power output and typical energy decrease with time after the 

■ initial few seconds of gamma-ray emission. This produces the X-ray afterglows, which have been detected up to 10 days after the 
burst (Frontera et al. 1998), and the optical ones, which have been detected up to a year after the burst (Fruchter et al. 1997, Bloom 
et al. 1998a, Castro-Tirado et al. 1998). 

The burst of May 8, 1997, was bright for a relatively long time and produced emission from gamma rays to radio. This enabled a 
detailed analysis of the expected spectral features of a synchrotron spectrum, confirming in great detail that we are indeed seeing syn- 
chrotron emission, and that the dynamical evolution of the expanding blast wave agrees with predictions if the blast wave dynamics 
are adiabatic (Galama et al. 1998a,b). In principle, one can derive the blast wave properties from the observed synchrotron spectral 
features. The problem is that the characteristic synchrotron frequencies and fluxes are taken from simple dimensional analysis in the 
published literature, so they are not suitable for detailed data analysis. Since there are now enough data on the afterglows of a few 
GRBs to derive their physical properties, we amend this situation in Sect. ||, correcting the coefficients in the equations for the break 
frequencies by up to a factor 10. We then use our theoretical results to infer the physical properties of the afterglows of GRB 970508 
(Sect. |J) and attempt the same for GRB 971214 (Sect.pJ). We conclude with a summary of results and discuss some prospects for 
future improvements in observation and analysis (Sect. ph. 

2. RADIATION FROM AN ADIABATIC BLAST WAVE 
2.1. Blast wave dynamics 
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We rederive the equations for synchrotron emission from a blast wave, to clean up some imprecisions in previous versions. Since 
the dynamical evolution of the blast waves should be close to adiabatic after the first hour or so, we specialize to the case of 
dynamically adiabatic evolution. This means that the radius r and Lorentz factor 7 evolve with observer time as (Rees and Meszaros 
1992, Meszaros and Rees 1997, Waxman 1997a, Wijers, Meszaros, and Rees 1997) 

r(t) = r dec (t/t dec y/ 4 (1) 

7(0 = ^Adec)" 378 . (2) 

Here r\ = E/Mqc 2 is the ratio of initial energy in the explosion to the rest mass energy of the baryons (mass Mo) entrained in it, the 
deceleration radius r de c is the point where the energy in the hot, swept-up interstellar material equals that in the original explosion, 
and fdec is the observer time at which the deceleration radius is reached. Denoting the ambient particle number density as n, in units 
of cm" 3 , we have 

, de c = (^^2) =1.81xl0'W«)' /3 Cc m (3) 

fdec = ^ =3.35(£ 5 2/») 1/ V(f s, (4) 

with OT p the proton mass and c the speed of light, and we have normalized to typical values: £52 = £'/10 52 erg and 77300 = ?7/300. 
Strictly speaking, we have defined n here as n = p/m v , where p is the ambient rest mass density. Setting f d =t/l day we then have, 
for t > f dec , 

r(t) = 2.29 x 10 17 (£ 5 2/n) 1/4 ^ /4 cm (5) 

7(0 = 6.65(£ 5 2/«) 1/8 f d " 3/8 . (6) 

Note that neither 7 nor r depend on r\\ once the blast wave has entered its phase of self-similar deceleration, its initial conditions 
have been partly forgotten. The energy E denotes the initial blast wave energy; it and the ambient density do leave their marks. It 
should also be noted that these equations remain valid in an anisotropic blast wave, where the outflow is in cones of opening angle 
9 around some axis of symmetry, as long as its properties are uniform within the cone and the opening angle is greater than I/7 
(Rhoads 1998). We should then replace E by the equivalent energy per unit solid angle £ = E/Q,. To express this equivalence we 
shall write the normalization for this case as £52 = £(4ir/10 52 erg), so we can directly replace £52 in all equations with £52 to convert 
from the isotropic to the anisotropic case. 

Before we can calculate the synchrotron emission from the blast wave, we have to compute the energies in electrons and magnetic 
field (or rather, summarize our ignorance in a few parameters). Firstly, we assume that electrons are accelerated to a power-law 
distribution of Lorentz factors, N(j e ) oc ■y^ p , with some minimum Lorentz factor 7 m . We are ignorant of what p should be, but it can 
in practice be determined from the data. The total energy in the electrons is parameterized by the ratio, e e , of energy in electrons to 
energy in nucleons. This is often called the electron energy fraction, but that term is only appropriate in the limit of small e e . The 
post-shock nucleon thermal energy is -fin p c 2 , and the ratio of nucleon to electron number densities is the same as the pre-shock value, 
which we can parameterize as 2/(1 +X), where X is the usual hydrogen mass fraction. In terms of these we have 



_ n e (E e ) _ 1+X m e p-1 7„ 
«7m p c 2 2 m p p-2 7 



£e - —3=-^ — — — (7) 



2 m„ p-2 

7m = T~T? — 7 £e7- (8) 

1 +X m e p-1 

The strength of the magnetic field in the comoving frame is parameterized by setting the field energy density, B' 2 /%-k, to a constant 
fraction, e B , of the post-shock nucleon energy density e' = 4--f 2 nm p c 2 . (Primed quantities are measured in the rest frame of the 
shocked, swept-up material; others are measured in the frame of an observer outside the blast wave at rest relative to the explosion 
center.) Consequently, 



B' = 7cy / 327rnm p e B 

= 2.58ey 2 £ 52 1/8 « 3 / 8 f d - 3/8 G. (9) 

From the above relations, we can express the evolution of the synchrotron spectrum from the blast wave in terms of observable quan- 
tities and six unknown parameters: £52, n,X,p, e e , and e B - But first we need to relate the synchrotron spectrum to these parameters. 



2.2. Synchrotron radiation 

We now derive the correct synchrotron frequencies and fluxes. These are strictly valid only for a uniform medium moving with 
a constant Lorentz factor. Real blast waves decelerate, of course, and have a more complicated structure behind the shock. The 
blast wave deceleration means that surfaces of constant arrival time are no longer the ideal ellipses expected for a constant speed 
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of the blast wave (Rees 1966) and at a given time we see contributions from gas with different Lorentz factors. This effect has 
been discussed thoroughly (Waxman 1997b, Panaitescu & Meszaros 1998, Sari 1998). The uniformity behind the shock is also a 
simplification: in reality the Lorentz factor varies from just behind the shock to the contact discontinuity. The density and other 
parameters vary accordingly (Blandford & McKee 1976). This effect has not yet been treated; it is expected to be comparable in 
importance to deceleration. Since both these effects are rather less important than our corrections to the synchrotron frequencies, we 
shall neglect both rather than attempt to apply only one of them. However, our improved treatment of the synchrotron emission is 
purely local and can be incorporated into any formalism that accounts for the varying local properties of the shocked medium at a 
fixed observer time. 

We assume that the electron population in any local volume has an isotropic distribution of angles relative to the magnetic field, 
and that the magnetic field is sufficiently tangled that we may average the emission properties assuming a random mix of orientation 
angles between the field and our line of sight. The radiated power per electron per unit frequency, integrated over emission angles is 

P'(V/V(7e), a) = ^ B F ( -. — \ erg cm -2 s _1 Hz -1 electron -1 , (10) 



m e c 2 \v±sma 

where F is the standard synchrotron function (e.g., Rybicki & Lightman 1979), and e and m e are the electron charge and mass, a is 
the angle between the electron velocity and the magnetic field, and 

37 e 2 <?fl' 

^(7e)=-p • (ID 

47rm e c 

(i.e., the traditional characteristic synchrotron frequency equals uxsina in our notation.) Next, we define the isotropic synchrotron 
function Fj S0 by averaging over an isotropic distribution of a. Setting x±("f e ) = u/uxije), we get 

. V3e 3 B' _, _, _, 

^isoC^x) = j- Fiso(xx) erg cm z s 'Hz 'electron (12) 

TO e C 

Fi so ( x ±) = / da sin aF(x±/ since) (13) 
Jo 

(We have made use of the symmetry of sin a to absorb a factor 1/2 into confining the integral to the first quadrant. The apparent sin- 
gularity at a = poses no problems because F decreases exponentially for large values of the argument.) Note that most calculations 
of blast wave spectra assume that the spectrum peaks at frequency j e 2 eB' / m e c. Due to the neglect of the factor 3/4-ir and the fact that 
F (x) peaks at x = 0.28587 and Fi so (x) at x = 0.22940, this estimate leads to quite erroneous inferences about blast wave properties. 

Finally, we must average the emission over a distribution of electron energies. We assume a simple power-law probability distri- 
bution of electrons between extreme values 7 m and 7 t : 

/(7e) = — ( — ) ' 7m < 7e < 7t (14) 



i-ht/j m y- p 

Now let x = u/uxQym)- Then the average power per electron becomes 



7m V7. 



. Vle 3 B' _ 7 _, _, _, 

Ppl( x ) = 5- Fpl(i) erg cm s Hz electron (16) 

m e c 

fo _d f e± 
Fpl(x) = —x 2 duu 2 F iso (u), (17) 

2 A7m77t 2 

in which we have transformed integration variable from 7 e to u = x~f m 2 /"f e 2 . The last equation shows the familiar result that for 
1 Si x < 7t 2 /7m 2 the spectrum from a power law of electrons is itself a power law. Since this region is known to extend over many 
decades in GRB and afterglow spectra, we quote numerical results for the case j x 2* 7m, for which the quoted results are independent 
of 7 t . The most easily identified point in the spectrum is its dimensionless maximum, x p , and the dimensionless flux at this point, 
Fpl(x p ) = 4>p> their dependence on p is shown in Fig.[l] Both now depend on the electron energy slope p. This defines the first two 
numbers that we can measure in the spectrum: 

, 3x„ •y m 2 eB' 

K. - (.9) 

Ui e C 

The calculation of the break frequency v c that separates radiation from slowly and rapidly cooling electrons (Sari, Piran, and 
Narayan 1998) is somewhat more difficult, because the cooling rate depends on both j e and pitch angle a. However, since the 
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cooling and the emission are both dominated by a = tt/2, we may estimate the break as the peak of F(x) for the value of j e where 
the cooling time for electrons with a = tt/2 equals the expansion time, t : 

6irm e c 

7c = (20) 
3 'v 2 eB' 

u e ' = 0.286^-^^- (21) 
4tt m e c 

In order to transform frequency and power from the rest frame of the emitting material to our frame, we note that the emission is 
isotropic in the rest frame by assumption. It is then trivial to compute the angle-average Doppler factors (see Rybicki & Lightman 
1979, Ch.4). For the received power, we find P = r y 2 (l+ 1 /3)P' ', which we shall simplify to P = 4j 2 P' /3 in keeping with the fact that 
our whole treatment is done in the ultrarelativistic limit, j3 — ► 1. Similarly, the intensity-weighted mean change in any frequency is 
v = 4"fi/' /3. Consequently, the appropriate mean of a power per unit frequency will transform as P v = ~{P' V . Of course, the spectrum 
also gets broadened, but that will not affect the locus of characteristic frequencies significantly. 

The synchrotron self-absorption frequency is usually set at the point where t v = 0.35. Using the co-moving width of the shock, 
Ar' = r/47, and the expression for the synchrotron absorption coefficient (Rybicki & Lightman 1979), we get 

3 IS 

^ = 47^/3 = 2.97 x 10 8 (^y iE^l(l+xf/ 5 n 3/5 e e - l e B l/5 S52 U5 (l+zr l Hz (22) 

where we have used equations ||-^ for the blast wave dynamics to express in terms of the unknowns we try to solve for, and added 
the correction for redshift, i.e. the equation in this form relates the observed frequency on Earth to the properties of the blast wave 
measured by a local observer at rest relative to the center of the explosion. Note that the self-absorption frequency in this simplest 
form is time-independent. We now also translate the other two frequencies into practical form: 

1^=471^73 = 5.73 x l0 1(, Xp (^) 2 e e 2 e B l/2 £52 1/2 (l+Xr 2 (l + z) 1/2 t d - y2 nz (23) 

v c =4*fv c '/3 = 1.12 x 10 12 e B -3 / 2 £ 52 ~ 1/ V 1 (l+zr 1/ V 1/2 Hz (24) 

Note the non-trivial redshift dependence of both, which stems from the fact that £<j is also measured on Earth and therefore redshifted. 
The observed flux at v m can be obtained by noting that our assumption of uniformity of the shocked material means that all swept-up 
electrons since the start contribute the same average power per unit frequency at v m (at any frequency, in fact), which is given by 
Eq. (|19|). Adding one factor of 7 to transform to the lab frame and accounting for the redshift, we have: 



N^PiJl+z) 



F "» = 4^ ■ (25) 



where N e is the total number of swept-up electrons, related to the blast wave parameters by N e = ^r 3 n(l +X)/2. The luminosity 
dista nce d epends on cosmological parameters, and for an £1 = 1, A = universe, which we shall adopt here, is given by di = 2c(l + 
Z-VT+z)/Ho. Consequently, 

h 2 

F Um = 1.15 , ™ ^M +X)£ S2 n x l 2 e B x l 2 mJy, (26) 



where h 70 = HJl km s" 1 Mpc" 1 . 

Equations [23| and g6| now are four independent relations between the four parameters of interest £52, n, e e , and e#. This 
means we can solve for all parameters of interest if we have measured all three break frequencies (not necessarily at the same time) 
and the peak flux of the afterglow. In addition this requires us to know the redshift of the burst, the electron index p, and the 
composition parameter, X, of the ambient medium. Note that multiple measurements of the same break at different times serve to 
test the model assumptions, but do not provide extra constraints on the parameters, since validity of the model implies that any of the 

four key equations is satisfied for all time if it is satisfied once. We therefore define the constants C a = u & /v &Jf , C m = iWdm / v m *, 
1 li 

C c = vj&c /v c *> an d Cp = Fi/ m /Fi/ m *. Here starred symbols denote the numerical coefficients in each of the four equations, and 
times denote the time at which the quantity in question was measured. Rearranging the four equations then yields 

3 

C-ic m -^C^C F hff p hp-^(^y-(\+X)-\l + zr- f Vi / ^~ 1 ) (27) 



"52 



5 11 1 1 -» ,1 1 / \-| 1 / \/\ +Z- 1 \ 

e e = C^C m --C c *C F ~x p "r P ( -^{^z) -(l+X)(l+zT-{ V j (28) 

e B = C-ic m -ic c ^C^xlf p hp-rr-^y 2 (l+X)(l+z)-i ( ^^ ) ( 29 ) 

_3 

n = C^cJ&c£c r -*xft<l>\(p-\Tl(*^~ i Q+^ (30) 

The last factor in each of these stems from the specific cosmological model adopted, and has entered the solution only via Eq. (^5]). 
To generalize to any cosmology, all that is needed is to replace (yT+z- l)/hjo in the above equations by (a! L /8.57Gpc)/Vl +z. 



5 



3. OBSERVED AND INFERRED PARAMETERS OF GRB 970508 

GRB 970508 was a moderately bright 7-ray burst (Costa et al. 1997, Kouveliotou et al. 1997). It was detected on May 8.904 UT 
with the Gamma-Ray Burst Monitor (GRBM; Frontera et al. 1991), and with the Wide Field Cameras (WFCs; Jager et al. 1995) 
on board the Italian-Dutch X-ray observatory BeppoSAX (Piro, Scarsi, & Butler 1995). Optical observations of the WFC error box 
(Heise et al. 1997), made on May 9 and 10, revealed a variable object at RA = 06 h 53 m 49?2, Dec = +79°16'19"(J2000), which 
showed an increase by ~1 mag in the V band (Bond 1997). BeppoSAX Narrow Field Instrument observations revealed an X-ray 
transient (Piro et al. 1997) whose position is consistent with that of the optical variable, and Frail et al. (1997) found the first GRB 
radio afterglow for GRB 970508; the radio source position coincides with that of the optical source (Bond 1997). 

The spectrum of the optical variable showed absorption lines at redshifts 0.77 and 0.835, indicating that 0.835 is the minimum 
redshift of the afterglow (Metzger et al. 1997). Subsequently, an [Oil] emission line with z = 0.835 was also found in the host's 
spectrum (Bloom et al. 1998a), which is often associated with star forming regions in galaxies. A faint underlying galaxy or star 
forming region is inferred to indeed exist from a levelling off of the light curve after 6-1 1 months (Bloom et al. 1998a, Castro-Tirado 
et al. 1998). It must be very compact, since the HST limits on an extended object underlying the GRB are fainter than the magnitude 
inferred from the light curve (Fruchter 1998). It is therefore almost certain that the compact nebula is the source of the [O II] line, and 
therefore also of the majority of the absorption lines. Given its compactness, a chance location of the burst far behind it is unlikely, 
and we shall assume that the burst occurred in this nebula, i.e. its redshift is 0.835. 

From the light curve behavior and broad-band spectrum (Fig. ||) of GRB 970508, Galama et al. (1998a,b) deduced the other 
properties of the burst required to calculate the physical parameters of the afterglow. We summarize them here: at t = 12.1 d after 
trigger, the break frequencies are f a = 2.5 x 10 9 Hz, v m = 8.6 x 10 10 Hz, and c c =1.6x 10 14 Hz. The peak flux is Fp m = 1 .7 mJy and 
the electron index p = 2.2. After the first 500 s electrons no longer cooled efficiently and the afterglow must evolve adiabatically. 
We shall set the cosmological parameters to be 51 = 1, A = 0, Hq = 70 km s" 1 Mpc" 1 . As noted above, they only enter the solution via 
the luminosity distance, and alternatives can therefore be incorporated easily via the substitution given below Eq. ([30|). Finally, we 
adopt X = 0.7 for the composition of the ambient medium. There are no reasons in the model to assume the ambient medium would 
not have normal cosmic abundance. While the metallicity Z is a strong function of redshift, X is hardly redshift-dependent, since the 
balance between H and He in cosmic matter has not been changed very much by nucleosynthesis. Using further that X2.2 = 0.580, 
02.2 = 0.61 1, we find 

£ 52 = 3.5 n = 0.030 

e e = 0.12 e B = 0.089. (31) 

We do note once more our deliberate use of £52, the energy per unit solid angle scaled to that of an isotropic explosion of 10 52 erg, 
in stead of the total energy: £"52 is truly constrained by the data, whereas the total energy requires us to know the as yet poorly 
constrained beaming of bursts. The recent findings by Fruchter et al. (1999) suggest there might be a break in the late light curve 
(100-200 days after the burst), at which time the Lorentz factor is 2 or less. If this break is due to beaming, it would be very modest 
beaming, and the total energy would be only a factor of a few less than the isotropic estimate. Our value of £52 does clearly rule 
out the very high energy estimates by Brainerd (1998) from the radio data alone. We have demonstrated for the first time that the 
electron and magnetic field energy densities are indeed close to but somewhat below equipartition. The ambient density is on the low 
side of normal for a disk of a galaxy but definitely higher than expected for a halo, lending further support to the notion that bursts 
occur in gas-rich environments. As an aside, we note that switching the values of v m and u c , which is allowed by the shape of the 
spectrum at 12.1 d, does not give a sensible solution (e.g., e e = 20.) This confirms the choice of Galama et al. (1998b), who noted 
that this solution was not compatible with the temporal evolution of the afterglow. 

The gamma-ray fluence of GRB 970508 was measured with BATSE to be (3.1 ±0.2) x 10" 6 erg cm -2 . Using z = 0.835 and h 70 = 1, 
this implies E^-y = 0.63. In other words, the energy emitted in gamma rays is 18% of the total blast wave energy (per unit solid 
angle in our direction). According to Galama et al. (1998b), the afterglow was cooling efficiently until 500 s after trigger; this means 
that during the gamma-ray phase all the energy given to electrons would be radiated away quickly, and mostly in gamma rays. If 
this phase is not too long, the energy radiated in gamma rays should be £s2 7 = e ej £ 52;, where e e7 is the value of e e during the early, 
gamma-ray emitting phase and £52; is the initial value of £52. Since the subsequent phase will be adiabatic, the blast wave energy 
measured at late times should be £52 = (1 - ee 7 )£s2i- Eliminating the initial energy, we conclude that 

C *7 _^52 7 (32) 



l-£e 7 f 52 

Therefore the measured ratio of gamma-ray fluence to late time blast wave energy implies that e e7 = 0.15, or slightly greater if some 
of the initial energy output is at E < 20keV. Compared with e e = 0.12 at late times, this demonstrates the near-constancy of the 
fraction of energy that is given to the electrons. Since the inferences about the initial gamma-ray fluence are independent of the 
whole machinery on blast wave synchrotron emission in the previous section, we may view this agreement as modest evidence that 
the coefficients derived there are close to correct, despite our simplification of the dynamics. 

It is also interesting to compare the properties at late times with those derived from radio observations. The scintillation size after 1 
month is about 10 17 cm (Frail et al. 1997), whereas our formulae give a transverse diameter of 5 x 10 17 cm; given the statistical nature 
of the scintillation size and our neglect of the gradients in properties in the transverse direction, to which this particular measurement 
is of course sensitive, this is not too bad. The Lorentz factor at this time is 3.4, so the evolution is still just in the ultrarelativistic 
regime. The field at this time is B' = 0.07 G. Katz and Piran (1997) estimated a size of the afterglow of GRB 970508 from a crude 
measurement of the self-absorption frequency. They found a size of 10 17 cm, and assuming an ambient density of 1 cm" 3 they found 
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that the Lorentz factor had already decreased to 2, and that most of the energy of the blast wave had been lost, i.e. it had evolved with 
radiative dynamics. The numbers we derive from our full solution after 1 week are: 7 = 5.8, transverse diameter = 2 x 10 17 cm. This 
means the blast wave is still quite relativistic, and with our low ambient density there is no need for radiative evolution. 

4. PROPERTIES OF GRB 971214 

This burst occurred on 1997 December 14.9727 UT. With a fiuence of 1.1 x 10~ 5 erg cm" 2 it is a moderately bright burst (Kippen 
et al. 1997). After localization by the BeppoSAX Wide Field Camera in X rays (Heise et al. 1997), the optical afterglow of this burst 
was found by Halpern et al. (1998). It shows evidence of strong reddening (Halpern et al. 1998, Ramaprakash et al. 1998). Once the 
afterglow had faded, a host galaxy became visible underneath it and its redshift was measured to be 3.42 (Kulkarni et al. 1998). 

One definite break was observed in the spectrum, and another possible one. The definite break (henceforth 'optical break') was 
found by Ramaprakash et al. (1998); they noted a break in the spectrum of the afterglow at 3 x 10 14 Hz, 0.58 days after trigger, in 
the extinction-corrected VRIJK spectrum. Another possible break (henceforth TR break') was found by Gorosabel et al. (1998) in K 
band data 3-5 hours after trigger. In Fig. [| we show the K band light curve. The data are not strongly inconsistent with a pure power 
law fit (x 2 /dof=2.0), but are suggestive of a break passing through K after about 5 hours. The physical interpretation of the afterglow 
depends rather strongly on whether the optical break is the peak frequency v m or the cooling frequency v c , so we shall discuss these 
two cases with their implications and problems in turn. 

4. 1 . The optical break as v m 

Ramaprakash et al. (1998) interpreted the optical break as the peak frequency, v m . A complication with the data is that the 
spectral slope is much too steep to correspond to any simple fireball model, which can be interpreted as due to reddening within 
the host galaxy (Ramaprakash et al. 1998, Halpern et al. 1998). Since reddening scales approximately as 1/wavelength, it cannot be 
determined without knowing what the true slope of the spectrum is. Assuming a blast wave model, one can predict this slope from the 
temporal decay rate of the flux and an interpretation of what break is seen. Following Ramaprakash et al., we now assume an adiabatic 
blast wave and interpret the break as u m . Then the flux above the break should depend on frequency and time as F = Fov~^t~ s , where 
(3 = 28/l>. There are several determinations of S, for the VRI fluxes: 1. 2 ±0.2 by Kulkarni et al. (1998)and 1.4±0.2 by Halpern et al. 
(1998). We shall adopt the value 1.3 ±0.2 in this paper, so (3 = 0.87 ±0.13 for the case under consideration. The left hand panel of 
Fig. |] shows the resulting de-reddened spectrum at t = 0.52 d (Note that Ramaprakash et al. construct the spectrum at a slightly later 
time, 0.58 d). We find that v m = 4x 10 14 Hz and F m = 30^Jy. (As an aside, we note that extinction fits with a mean galactic extinction 
curve are worse, since the redshifted 2200 A bump falls within VRI.) The amount of extinction, 0.43 mag at a rest frame wavelength 
of 5500 A, is very modest and does not imply a special location of the burst within the galaxy. A strong point of this fit is that the 
X-ray flux measured at the same time, which was not included in the fit, agrees nicely with it. The reported non-detections in the 
radio at levels of 10-50 /iJy could be inconsistent with the peak flux: the reddening-corrected F m is 30 /iJy, so we may have to invoke 
self-absorption to suppress the flux at 8.46 GHz. A weak point is the fact that the flux at K, which is below the peak frequency at 
0.58 days in this model, would have to rise with time as t 1 ' 2 as the peak approaches. But the early K data by Gorosabel et al. (1998) 
clearly indicate a significant decline of the flux in K from 0.21 to 0.58 days. In fact, the early K flux even exceeds the supposed peak 
flux in the spectrum at 0.52 d, which means the peak flux had to be declining as well. This requires a non-standard blast wave model 
(e.g. a beamed one or a non-adiabatic one) and is thus inconsistent with the blast wave model used in this interpretation. Nonetheless, 
we shall briefly explore the physical implications of this model, since the interpretation of the optical break as a cooling break is not 
free of problems either. 

For the simple adiabatic model used by Ramaprakash et al., we get the blast wave energy from Eq. (Q): 

£^=21^(^) A (^)\ (33) 
l+zVO.2/ V0.1/ 

The coefficient is 60 times larger than in Eq.3 of Ramaprakash et al., almost solely due to our more accurate calculation of the peak 
frequency. We can use Eq. (E6) to derive an independent estimate of the blast wave energy from the peak flux of 30/^Jy: 

£ 52 = 0.09»- I / 2 (^-)" I/2 (34) 

This value is difficult to reconcile with the energy estimate from v m , unless we push the equipartition fractions very close to unity 
and/or adopt a very low ambient density. 

4.2. The optical break as v c 

In order to accomodate the decline in the early K band flux, we now assume that v m was already well below K band at 0.52 d, 
and therefore the optical break is v c . Then the spectral slope in VRI should be related to the temporal decay at those frequencies by 
f3 = (26+ l)/3, which for S = 1.3 ± 0.2 implies f3= 1.2 ± 0.13 and p = 2.4. This model is explored in the righthand panel of Fig. El 
While it solves the K decline problem, it is a worse fit to the K flux at this time, and does not do very well in predicting the X-ray 
flux, which is 23<r above the extrapolated spectrum. Also, if we then say that the IR break is real, we have a peak flux of 60 /iJy. 
This value is greater than the 8.46 GHz flux limits obtained 0.8-20 days after the burst (Kulkarni et al. 1998). This would require 
that the self-absorption frequency in this afterglow exceeds 10 GHz. Alternatively, the peak frequency at 20 days could be at least 
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a factor 200 above 8.46 GHz, so that an extrapolation from the peak flux to the radio using an optically thin synchrotron spectrum 
(Fi, oc i/ 1 / 3 ) falls below 10 /iJy. Since we know in this case that v m was 1 .4 x 10 14 Hz after 0.21 days, we can use v m oc t~ 3 / 2 to find 
that v m = 150 GHz at 20 days, too low to be compatible with the radio upper limits; we conclude that we must require v & <: 10 GHz, 
as in the other interpretation of the break, to suppress the radio flux. 

Using v m and F m from the IR peak at t = 5 h= 0.21 d, we can again use equations ^ and ^6| to get two expressions for the energy 
in terms of the other unknowns: 



e B y 1 



= 3 -°(o^r(o59j (35) 

Ss2 = l.l(^V I/2 (^r 1/2 (36) 
V 0.030/ V 0.089/ 

Here we have scaled the unknowns to the values found for GRB 970508. In this case, the two independent energy estimates are quite 
compatible. 

Now that we identified all but the self-absorption frequency in the afterglow of GRB 971214, we may use equations p7^-BQ to get 
all the parameters of the burst, leaving their dependence on the unknown v a explicit. The cooling frequency is v c = 4 x lC^Hz at 
0.52 days. It follows that 

£52 = 0.46^ » = 0.60^ 6 Z 
e e = 0.26^ z e B = 0.027 v~^. (37) 

Since we require f a > 10 GHz in order to satisfy the radio limits, we get £52 < 0.07, n > 9 x 10 4 , e e > 1 .8, and e B < 8 x 10~ 5 . These 
values are rather different from those of GRB 970508, and e e > 1.8 is implausible (but not impossible). With such a low energy and 
high density, the GRB would become non-relativistic (and start a faster decline) within about 4 days. This means that the late time 
radio flux may be suppressed without invoking v. d > 10 GHz, easing the constraints on the parameters somewhat. But with v. i =\ GHz 
we get values for which the GRB remains relativistic for 400 days; this conflicts with the radio limits. We conclude that 24 probably 
has to exceed 5 GHz for a consistent solution, implying a fairly low energy and high ambient density for GRB 971214. With a high 
emitted gamma-ray energy per unit solid angle, 8si 1 = 30, this means that the ratio of emitted gamma-ray energy to remaining energy 
in the blast wave is uncomfortably large, probably in excess of 100. Either very efficient radiation in the GRB phase of the burst is 
needed to achieve this, or a strong difference in the amount of beaming between the gamma-ray and afterglow emission. Beaming 
of GRB 971214 to a degree similar to that observed in GRB 990123 (Fruchter et al. 1999, Kulkarni et al. 1999) would allow a more 
standard solution for the parameters: all the trouble arises from the fact that the radio limits force a high v. d . If the burst was beamed 
so that it started declining more sharply after a day or so, then the radio flux at later times could be naturally low, and we would be 
allowed to choose v a ~ 1 GHz, which leads to parameter values for GRB 971214 that are similar to those of 970508. It would even 
open the possibility that the temporal decay is somewhat contaminated by the steeper decay after the beaming break, so we should 
have used a slightly shallower temporal and spectral slope in fitting the spectrum of GRB 971214 at 0.52 d with a cooling break. Then 
both the X-ray and K flux would agree much better with the fit. However, this now is a little too much speculation for the amount of 
data available, and it may well be wrong to try to bring the parameters of GRB 971214 into agreement with those of GRB 970508. 
We already know that afterglows are diverse. For example, GRB 980703 had a cooling break in X rays after a day rather than in 
the optical (Vreeswijk et al. 1999, Bloom et al. 1998). We simply need more well-measured bursts in order to establish the allowed 
ranges of parameters like e e and e B . 

In summary, neither the identification of the optical break as v m nor as v z is without problems, so we should take any derived 
parameters for this burst with a grain of salt. But the v m interpretation is in our view the least tenable, since it predicts that the flux 
at K should have risen from the start until about a day after the burst, whereas the data clearly show a declining K flux. The problem 
of the v c interpretation is that the X-ray flux at 0.52 days is about 2.3er higher than expected, and the K band flux lower, which is 
somewhat uncomfortable but perhaps tolerable. 



5. CONCLUSION 

We have calculated the synchrotron spectra from the blast waves causing GRB afterglows and derive improved expressions for 
the relations between measured break frequencies and the intrinsic properties of the blast wave. These allow us to relate the blast 
wave properties to observable quantities more accurately. We correct the expression for the blast wave energy by almost two orders 
of magnitude. Our expressions are exact for an undecelerated, uniform medium. Deceleration and radial structure of the shock are 
expected to change the expressions for the final parameters by another factor few, much less than the corrections found here but 
still of interest. Combined with the uncertainties in the measured values of the spectral breaks and fluxes, this means that the blast 
wave parameters derived here are still uncertain by an order of magnitude (see the solution by Granot, Piran, and Sari 1998 as an 
illustration of the possible differences). 

There is enough data on GRB 970508 to compute all intrinsic parameters of the blast wave. The energy in the blast wave is 
3 x 10 52 erg/47r sr. The ambient density into which the blast wave expands is 0.03 cm" 3 , on the low side for a disk of a galaxy. The 
fraction of post-shock energy that goes into electrons is 12%, and that into magnetic field, 9%. We also estimate the fraction of 
energy transferred to electrons during the gamma-ray phase, and find this to be 15%. The agreement with the later blast wave value 
suggests that the fraction of energy given to electrons is constant from 10 s to 10 6 s after the trigger. 
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For GRB 971214 there is ambiguity in the interpretation of the break seen in the optical half a day after the burst. We argue that 
the break is most likely to be the cooling break, but the argument is not watertight. Assuming it is the cooling break, we still lack 
the self-absorption frequency, but radio limits constrain this to be in excess of about 5 GHz. The limits on parameters that follow 
from this indicate that the afterglow properties of GRB 971214 are different from those of GRB 970508. GRB 971214 must either 
have had more narrow beaming in gamma rays than in optical or have radiated its initial energy with more than 99% efficiency in the 
gamma-ray phase, according to the parameters we derive. Also its magnetic field was far below equipartition, e# < 10~ 4 . However, 
beaming or other additions could ease the constraints, and allow parameter values similar to those of GRB 970508, so the physical 
parameters of GRB 971214 are very poorly constrained. 

Our analysis emphasizes the importance of early measurements covering a wide range of wavelengths. The full identification of 
the cooling frequency v c in GRB 970508 hinged on abundant photometry, including colors, being available soon after the burst, since 
the break passed R after 1.5 days (Galama et al. 1998b). In H and K, the action lasted a week (Galama et al. 1998b), and this is the 
general trend: there is more time in IR, since all breaks pass later there. However, our revised coefficient for the peak frequency, u m , 
shows that the peak can only be caught in the IR within hours of the trigger (or much later in the radio). A case in point are the very 
early K' band measurements of GRB 971214 by Gorosabel et al. (1998), which provide an invaluable constraint on this afterglow as 
they may have caught the passage of v m through K' . Therefore, we encourage first and foremost early long-wavelength coverage, 
including searches for afterglows in IR, as a method of effectively constraining afterglow parameters. Two of the three crucial break 
frequencies in an afterglow can pass the optical and IR within hours and days, respectively. There is no time to first search and 
only then attempt broad coverage. Instantaneous alerts from HETE2 and SWIFT will therefore greatly advance our understanding 
of afterglow physics. For HETE2, and to some extent for SWIFT, this will require an enormous amount of work from a network of 
ground-based observatories with good coverage in longitude and latitude, so that always at least one observatory is well-placed for 
immediate response. 



R. A.M. J. Wijers was supported for part of the work by a Royal Society URF grant. T.J. Galama is supported through a grant from 
NFRA under contract 781.76.01 1. We thank Jules Halpern for alerting us to the X-ray data on GRB 971214 and pointing out that our 
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P 

Fig. 1.— The dimensionless location x p (solid) and dimensionless peak flux <f> p (dashed) of a synchrotron spectrum from a power 
law of electrons, as a function of the power law index, p of the electron energy distribution. 




Fig. 2.— The X-ray to radio spectrum of GRB 970508 on May 21.0 UT (12.1 days after the event) from Galama et al. (1998b). 
Indicated are the inferred values of the break frequencies v- d , v m and v c for May 21.0 UT. 
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Fig. 3. — The K flux of GRB 971214 as a function of time. Extrapolations to early times of the Keck measurement after 14 hours 
(Kulkarni et al. 1998) are shown for the same three values of 8 used in Fig.Q. Other data are by Gorosabel et al. (1998) and Garcia et 
al. (1998; the upper limit). The assumption here is that K lies above v m and below is c (the situation in the righthand panel of Fig. Q). 
Therefore, the actual slope of each curve is 0.25 less than the value of 8 (by which it is labelled) because K lies below the cooling 
break, and <5 is measured in VRI, above the cooling break, where the temporal slope is steeper by 1/4 than below it. 



Broadband spectrum of GRB971314 at Dec lo.DO UT Broadband spectrum of GRB971214 at Dec lb. 50 UT 
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Fig. 4. — The near-infrared/optical to X-ray spectral flux distribution of GRB 971214 on Dec 15.50 UT (the epoch of the / band 
measurement). Open symbols indicate the measured values (from Ramaprakash et al. 1998 and Heise et al. 1999) extrapolated to 
Dec 15.50 UT. Note that the error on the J band data is much larger than used by Ramaprakash et al. (1998), in agreement with 
the original report (Tanvir et al. 1997). Filled symbols indicate the dereddened data. The dotted lines indicate the 1-a errors on the 
spectral slope as derived from the temporal slope, (left) Result of assuming that the break is the synchrotron peak. The spectrum 
below the peak follows the low-frequency tail of the synchrotron spectrum, where F v oc u 1 ' 3 . (right) Result of assuming that the 
break is the cooling break. The spectral slope changes by 0.5 across the break. The extinction, Ay, derived from the fit, corresponds 
to the rest frame V band. 



